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SUMMARY 


The purpose of this research has been to design and develop a flexible 
source of spinning modes which is capable of generating independent spinning 
waves of controlled complexity and spin speed without the introduction of 
broad band elements. These features were accomplished through the use of 
eight commercial loudspeakers located in an equally spaced circular array 
with diameter of 11 inches and properly phased so that the system could 
generate a spinning wave. The constructed apparatus was tested in an anechoic 
environment and found capable of generating a plane, one and two lobed spin- 
ning wave of high quality with a sound pressure level of 120 dB and at frequen- 
cies ranging from 1500 to 2500 Hz at a distance of V in the far field. The 
wave speeds investigated varied from 8000 to 18000 rad/sec which represent 
supersonic peripheral speeds. 


INTRODUCTION 


Today it is the accepted belief that the periodic redistribution of 
pressure over the profile of fan or compressor rotor blades as they intersect 
the wakes of upstream objects can give rise to the phenomenon of propagating 
spiral waves. In the past, research on developing an effective and economical 
design for attenuating these spiral waves by absorption and cut-off techniques 
has been frustrated by the presence of a severely inhomogeneous turbulent 
velocity and temperature spectra through which the wave must propagate. In 
particular, those studies directed toward relating pertinent properties of 
duct lining configurations operating in a jet engine service environment are 
further complicated by the presence of a multitude of these spiral waves along 
with a high level of background broadband noise. In general, the duct geometry 
is highly complex and sound levels are sufficiently high to warrant the need 
for using the nonlinear wave equation to describe propagation. 

Thus the possibility for a meaningful parameter variation that could 
lead to a more fundamental understanding of the characteristics of propagating 
spiral waves in actual turbomachinery applications is limited. Hence, there 
has arisen a genuine need for the development of a flexible and predictable 
source of spinning waves so that the pertinent parameters characterizing 
their propagation can be more easily related to absorption properties of 
various configurations of different lining materials. In addition, the propa 
gating characteristics of spinning waves in ducts of varying cross sections 
is an area of considerable interest and could be studied more effectively. 

The objective of this research was the design and development of such a device. 
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THE ANALYSIS 


WAVEFORMS TO BE SYNTHESIZED 


As a first approximation we shall assume that the propagating waveform 
can be described by linear acoustics. Therefore, away from the source region 
the complex wave structure of these spinning waves may be effectively synthe- 
sized by a linear combination of normal modes each of which satisfy the 
homogenous wave equation and prescribed boundary conditions. The selection 
of the boundary conditions ultimately fixes the form of the normal modes 
as well as the eigenvalues in the expansion. Thus selection of the boundary 
conditions determines the form of the proposed Spinning Mode Synthesizer. 

In order to select an economical and representative design, we chose 
for the boundary conditions a straight, hard-walled concentric cylindrical 
duct terminated anechoically . This permits the use of boundary conditions 
describing an infinite duct, or that of acoustic waves propagating in one 
direction only. The purpose for doing this is that the forcing function in 
the non-homogenous wave equation only serves to specify the numerical values 
of the Fourier components in the eigenfunction expansion. The composite 
form of these Fourier amplitudes will depend upon the particular engine 
design. To date, there has thus far been no adequate analytical treatment 
for evaluation of these coefficients, and until there is, one must rely upon 
experimental measurements. But, as will become clear later, it is always 
possible to adjust the amplitude and phase c£ these components to any desired 
level. Hence, we seek only to develop a device capable of generating indivi- 
dually the discrete set of eigenfunctions that could contain propagating 
acoustic energy in the duct of a rotor. We are now in position to investi- 
gate the form of the discrete set of eigenfunctions that are allowed subject 
to the above assumptions. 


ALLOWED COMPRESSOR MODES 


The first penetrating analysis on the form of the permitted normal modes 
and the criterion for their propagation or decay in terms of pertinent com- 
pressor variables was presented in Tyler and Sofrin’s [l] classic paper on 
’’Axial Flow Compressor Noise”. 1 Basically, they obtained their result by 
analyzing the acoustic field away from the source region. It was, therefore, 
possible to use the standard technique of separation of variables to the homo- 
genous wave equation with the previously mentioned boundary conditions. This 
permitted a solution for the form of the eigenfunctions that could be used 
in a Fourier synthesis of the sound field. Their result is presented as 
equation (1), and the coordinate system is defined in Figure 1. 


P(r,0,Z,t) = E Z 

n=0 m=-° 3 


[A cos (mO - cot + 
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where 


E° (k° r) 
my my 


J (k r) 
m my 


J ’ (abk , 
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Y m' (0bk my) (k ^ r) 


and 


J^, are cylindrical Bessel and Newman functions respectively 

k^, m> k z > w are eigenvalues of the (r, 6, Z, t) coordinates respectively 

A , <b are amplitude and phase coefficients of the mnth wave 

mn" r mn r r 

y radial mode index 

o hub-tip ratio (a/b) 

m circumferential mode index 

n harmonic of blade passage frequency 

where primes indicate the first derivative with respect to r. 

In order to determine the specific form of the eigenvalues, the boundary 
conditions must be applied. However, since only three of the eigenvalues are 
independent, Tyler and Sofrin chose the Z coordinate eigenvalue to be repre- 
sented in terms of the characteristic time and radial modes 


0 ). 2 


k = [(-) 


(k° ) 2 ] 1/z 
v my' J 


( 2 ) 


The radial boundary condition for a hard-walled duct requires that 


^a,b E 0 


and application of these boundary conditions force to satisfy the 
relation 


my 


J ' (bk° ) Y ' (0bk° ) - J ' 
m v my' m v my' m 


(abk a ) Y ' (bk a ) = 0 
my' m ^ my 7 


(3) 


In their analysis Tyler and Sofrin do not describe how one is to obtain the 
(0,t) eigenvalues. In order to satisfy the solution of the equation, the 
allowed order of circumferential modes can be shown to be given by equation (4). 
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(4) 


m = n B + h V 


The reader is referred to Appendix B, where a derivation is given for obtain- 
ing such a relation. Since the interaction frequency (nBfi)-l is available 
for forming a characteristic time it is reasonable to assume that oa = nBQ. 

In their paper, Tyler and Sofrin present an excellent discussion of why 
sound modes of increasing complexity, increased waviness, are subject to 
decay in the duct and why lower order modes can propagate without decay. The 
condition when — = is known as duct cut-off. This condition is reached 
by essentially considering whether the value of k z is real or imaginary in 
equation 1. Decay occurs when k z is imaginary, and it can be shown from 
equations 2 and 3 that this tends to be the case when m and y increase in 
value. Thus it appears that the most fundamentally important modes to investi- 
gate from the standpoint of unattenuated propagation would be that of lowest 
order. We have adopted this procedure in the design of the experimental 
apparatus by limiting the number of sound sources in the final array. 

It is worthwhile to remember that equation 1 only serves as a model 
for structuring the sound pressure distribution in actual turbomachinery 
applications. One must eventually take into account the inlet flow as well 
as the duct ? s finite length and termination geometry. Since precise speci- 
fication of the rotor field's acoustic source distribution is still not 
possible, the spectral energy level and distribution in wavenumber space for 
propagating duct modes remains difficult to relate to various aerodynamic 
variables. Therefore in the development of the synthesizer we would have to 
assure that all propagating modes contain an equal amount of acoustic energy. 
This restriction gives rise to the design criterion that the geometrical con- 
figuration of the proposed device should produce the flattest possible output 
spectrum in wavenumber space. 

Thus in an attempt to provide a greater flexibility for the synthesizer, 
the sound pressure field in a semi-infinite duct terminated at one end was 
computed for a source ring distribution as shown in Figure 2< The details 
of this calculation are given in Appendix C and the results of this calcu- 
lation (Equation C.6) is reproduced here as equation 5. As one can observe 


P(x,t) 
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from the form of equation 5, even though the dipole sound source distribution 

at z'=0 is selected to produce one particular 0-mode, all possible radial 

modes are also excited. However, each one is diminished by the Jm( k my r °3 

term appearing in equation 5. By using the tables from Jahnke and Emde [4] 

and the radial boundary condition J ' (k ) I = 0 to compute the k 's, the 

J m my ' r my 
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amplitude distribution of this term appears in Figures 3, 4, and 5 for • 
the source ring radius as a percentage of the duct radius. Since only one 
0-mode is developed at one time, depending upon the source distribution, 
each figure is thereby cross-plotted with the first few radial modes from 


0 - 3. 


As we have discussed previously, each (m,y) mode becomes propagational 
at a different source frequency depending upon whether k 2 = [k 2 -k m y 2 ] */ 2 
is real or imaginary. The frequency at which a mode becomes propagational 
is known as its cutoff frequency. Each cutoff frequency for the radial 
modes appearing in figures 3,4, and 5 is listed to show, for example, for a 
12"dia. duct the number of modes which can propagate depending on the source 
frequency selected. In general, however, what these curves indicate is that 
there is really no one optimum location for a source ring if one also desired 
control over the first few radial modes. Realistically though, one must choose 
a fixed source radius. The most efficient coupling to the duct apparently 
occurs between 10 to 20% and 90 to 100% of the duct radius. In our design we 
chose to locate our sources at 0.96 b. While this fixed radial position did 
not give us the greatest flexibility as compared to a variable ring radius, 
it does provide for the maximum number of (m,y) modes to be generated effi- 
ciently. By using a variable source ring radius and the criterion for cutoff, 
one could conceivably selectively generate and propagate only one (m,y) mode. 

The results of this idealized analysis also indicates where not to 
locate the source ring. For the first few (m,y) modes desired, care should 
be taken to avoid the zero crossings of the curves in figures 3,4, and 5 if 
one is to select a fixed source ring radius. Obviously, at these points an 
infinite source energy is required to propagate energy in the corresponding 
mode. 


At this point then we must demonstrate that a discrete set of acoustic 
sources can generate one discrete 0-mode when coupled to a duct, realizing 
that the proper radial spacing of the sound sources is adequately determined 
as in Appendix C. Therefore equation 6 represents the (0,t) 


P(0t) = A cos(m 0 - nBQt + <6 ) 

v mn v Y nur 


( 6 ) 


distribution of sound to be modeled from equation l,and indicates that there 
will be for each harmonic of blade passage frequency an infinite number of 
left and right running lobed waves traversing in a helical path up the inlet 
duct. The waviness of the allowed lobed waves is determined by the index m. 
Figure 6 is a visualization of the concept of the spiral mode as it propa- 
gates down (or up) the duct and figure 7 represents for a fixed time the 
sound pressure distribution in the (r,0) plane. 

Since one of the major design features of the spinning mode synthesizer 
is the ability to generate one lobed structure at a time corresponding to a 
single blade passage frequency, equation 6 may be reduced to 


P(0,t) = A^ sin(m 0 - nBQt) 


(7) 
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where 4> m h is dropped because it represents only the phase relation among 
waves. Equation (7) corresponds to the acoustic pressure of one wave 
located in a plane normal to the duct axis. It is this relation which we 
wish to simulate by the spinning mode synthesizer. 


ARRAY EQUATIONS 


The preliminary objective in the development of the spinning mode 
synthesizer in this research was to demonstrate the feasibility of satis- 
fying equation 7 in a free field test. Later developments were directed 
toward satisfying the conditions for a hard-walled duct. The easiest way 
to develop circumferential waves as defined in equation 7 would be to 
construct a device where sound sources are distributed circumferentially in 
a planar array as shown in figure g. In order to satisfy equation 7, the 
drivers are sequentially phased around the circumference of the array. With 
the use of the coordinate system shown in figure 8, it can be shown that a 
system of n piston sources of equal strength (see Morse, 1968 [6]) radiates 
sound to the free far field as given by equation 8, 

n 2J i (ka sin u.) sin(kR.-a)t + <f>.) 

P(x,M,t) - ACk.D) I taslnu . J i- (8) 

J =1 J J 


where A(k,D) source sound pressure amplitude 
a driver radius 

k acoustic wave number of the driven piston 

D driver separation 


R. , U. 
J J 


distance and solid angle respectively to radiation position 
from each source as defined in figure 8. 


When u is small, it is possible to choose a piston diameter of sufficient 
size to generate ample sound intensities but also so that 2Ji (ka sin u.)/ 
ka sin U.-l. Also in the far field small u. implies that R. is practically 
constant^ Therefore, besides the amplitude^being a functioA of frequency 
and relative separation between sound sources, the sound pressure distribution 
in the far field is primarily sensitive to 


n 

E [sin(kR. + <p. - cot] (9) 

j=i 3 3 


The analogy with 


equation 7 becomes complete when we adjust <J> . such that 
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sin(m0-nBftt) 
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Z 

j = l 


[sin kRj + ^2. (j-1) - cot] 


( 10 ) 


In other words, we should be able to generate spiral type waves when we 
provide phasing between drivers based on the simple relation 


360m 


(j-1) 


(ID 


when developing the (m,n)th wave. 

In order to demonstrate that it was at least possible to develop the 
desired waveforms using the piston driver as a model, equations 8 and 11 
were used to compute the free field sound pressure for various values of 
the parameters (m,f,i^,n,D). In every case the desired waveform agreed 
identically with the expected distribution of sound pressure for the 
selected values of the above parameters. Later we will discuss the criterion 
for their selection, but numerically the circumferential mode index m was 
varied from 0 to 4, the driver frequency f from 1000 to 2500 Hz in 500 Hz 
increments, and the solid angle ^ from 10° to 20° in 5° increments. The 
number of drivers and their linear spacing D was fixed at 8 and 4.25" res- 
pectively. 

Typical results of these calculations are shown in figures 9 through 
lg for f = 2000 Hz only. As is immediately obvious from the curves of 
nondimensionalized sound pressure versus circumferential angle 0, the distri- 
bution of pressure is as expected for modal shapes 0 through 4. Perhaps 
though the decrease in SPL with a corresponding decrease in ip is not as 
apparent. As will be demonstrated later, this occurs for all modes except 
m=0, since theoretically complete cancellation occurs at t|j= 0. This is a 
property of free-field operation only, for when coupled to a duct, the 
radial distribution is determined soley by the boundary conditions. 

The time history at any two points, 0 and 0+60, provide the opportunity 
for observation of the phase velocities of the generated waveforms. The 
curves in 3.4 through 18 :for SP vs t clearly indicate each mode's direction of 
spin and angular velocity. For the jet engine compressor, the duct mode spin 
velocity is given by nBQ/m. For our proposed synthesizer, the modal spin 
velocity is given by w/m as can be determined from these curves. If the base 
frequency is tuned to any multiple of the blade passage frequency so that 
(j = nBQ, then our model's wave speed corresponds exactly to that of an actual 
engine for the same modal pattern. As is shown in the figures, the spin 
direction is determined by which of the two curves is leading or lagging the 
other in as much as the spatial point of measurement is known. 

With regard to the feasibility of the spinning mode synthesizer, the 
proposed model performed satisfactorily. With appropriate phasing, high 
quality circumferential modes ranging from 0 to 4 can be generated under the 
idealized conditions as imposed by equation 8. Not only can the wave speed 
be controlled and determined, but also the wave's spin direction. Therefore, 
at least theoretically, the proposed model appears to be the flexible source 
of higher order duct modes for which we were originally in search. 
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THE EXPERIMENT 


EXPERIMENTAL APPARATUS 


In the last section, one sees that by using eight piston drivers with- 
out considering the effect of near field cancellation, one is able to 
generate the desired waveforms. In the design of the synthesizer, however, 
we had to provide for the most efficient arrangement of commercial loud- 
speakers which would minimize near field cancellation and at the same time 
behave as a piston type sound driver. 

With regard to the acoustic efficiency of the array, Pritchard [7] 
presents curves of the normalized mutual radiation impedance for two piston 
drivers as a function of relative separation. This is his figure 2. By 
selecting D = 6", a = 1/2" and varying f from 1 to 3 kHz, Pritchard*s data 
predicts that the mutual radiation impedance would be definitely within 
+0.2 for drivers in phase. While use of his data for drivers not in phase 
may be questioned, the curves still provide valuable guidelines to the 
design of a proposed system for minimizing near field cancellation. 

With the spatial arrangement of the sound drivers determined by Pritchard f s 
data, we found a rather inexpensive but extremely adequate loudspeaker with 
good frequency response and directionality characteristics similar to that 
of a piston driver. Figure 19 shows the selected loudspeaker and the agreement 
obtained between the two for loudspeaker frequencies of 1, 2, and 4 kHz. Since 
the device was to later be mounted to a duct, the solid angle ^ was limited 
to 0 to 20°. As can be seen for 2000 Hz , exact agreement is obtained between 
the piston driver and the selected loudspeaker directionality patterns for 
ip < 20°. One, therefore, can expect reasonable operation of the proposed design 
in a duct for it satisfies the basic conditions of the analytical assumptions. 

The experimental apparatus that was decided upon is shown schematically 
in figures 20 through 22. The overall assembly of figure 20 was mounted in 
Penn State* s Anechoic Chamber for all subsequent free field tests. The 
sound driver array plate, figure 21, was installed in a 4 ft. diameter 
wooden baffle to simulate radiation to one side. The rear of the array was 
enclosed with a 2 ft. diameter fiberglass filled container. The purpose of 
this being to load each sound driver uniformly. Since each driver *s output 
varied differently as a function of its input, it was found necessary to 
construct a probe tube microphone which could be rotated remotely to each 
sound driver. The output phase and sound pressure level of each driver could 
then easily be adjusted to desired values. 

Two microphone processing systems were used for measuring the generated 
sound field in parallel and perpendicular planes from the driver array. The 
vertical plane of figure 8 was measured by the system shown in figure 22. 

As can be seen a 1/4" condenser microphone is mounted so that it could be 
rotated through 360° and indexed from 1** to 17" radially and up to 4* - 6* 
away from the sound driver array. The horizontal plane was processed by using 
a remotely controlled boom microphone system. A 1" condenser microphone was 
used to record the directivity patterns of the entire array from radii of 2.5 
to 4 feet. 
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EXPERIMENTAL PROCEDURE 


The purpose of the experiment was, of course, to verify that one could 
selectively generate higher order 0 duct modes with the proposed design. In 
order to demonstrate this, measurements in planes parallel to the driver 
array were carried out according to the following schedule. Since very little 
difference exists between the radiated sound patterns for the ducted and un- 
ducted driver array when the solid angle is less than 20° (See figure 8), 
all of the following preliminary testing was conducted with the array mounted 
under free field conditions. The 1/4" microphone was positioned at 2.5, 3.0, 
3.5, and 4.0 ft. away from the driver plane and rotated 360° for microphone 
radii from 1" to 17" which corresponded to the solid angles of 5°, 10°, 15°, 
and 20° at each respective location. This was done for modes 0 through 3 
and for driver frequencies of 1500, 2000, and 2500 Hz . Figures 23 and 24 
characterize the electronic drive system used to supply each of the eight 
drivers in the array and also the electronics used in processing the data from 
the 1/4" and 1" microphones. 


EXPERIMENTAL RESULTS AND DISCUSSION 


We will first present data for the array* s directivity. In figures 

25 through 28 solid lines show the recorded directivity of the entire 
array for a driver frequency of 2000 Hz only for modes 0 through 3. The 
only observed change when the driver frequency increased was a decrease in 
overall beam width for all patterns. The dashed curves are generated from 
Tyler and Sofrin's equation 7.3.2 which describes radiation from a circular 
duct with the additional possibility of radial modes. It is particularly 
interesting to note the remarkable agreement between the array* s directivity 
and that predicted by 7.3.2. The agreement appears best with the radial 
mode y = 1. This corresponds to physical reality, since if we recall figure 

26 from Tyler and Sofrin, we observe that if we phase our drivers for m =1, 
the mechanical location of the drivers dictates approximately a near field 
pressure distribution of y = 1 as shown in the figure. It is not surprising 
then that with a single ring of drivers the near field pressure pattern as 
arranged in our design would always simulate a (m,l) mode of an open face 
circular duct. 

Because of the very good agreement exhibited in figures 25 to 28, one should 
be able by using various number of driver rings to simulate accurately any 
(m,y) mode radiated from a fan or compressor inlet duct. The utility in doing 
this may be realized if one seeks to determine the original near field pressure 
pattern of a frequency spectrum from the noise field of a full scale engine. 

Such a spectrum appears in figure 26 from directivity measurements carried 
out by Morfey and Dawson [8] on a full scale turbojet engine. The data used 
corresponds to their engine B and the directivity measurements of figure 2a 
where in our case the S/s A parameter equals 0.393. Agreement appears good up 
to 45° off axis, where beyond this the turbo jet*s broadband spectrum becomes 
dominant . 
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Perhaps more convincing results are presented in figures 29 through 
35. These curves represent the measured RMS sound pressure in a plan para- 
llel but *4' distant from the sound driver array. For the sake of brevity, 
only the data taken at this distance for modes 0 through 3 at a driver 
frequency of 2000 Hz are presented. There are actually two sets of curves. 
Figures 25 through 28 delineate the measured phase angle between positions 
in the plane for microphone radii of 7.625”, 12,345”, and 17.275”. The 
second set of curves, figures 29 to 35, represent for any arbitrary fixed 
time the resulting spatial waveform. 

An obvious bounding condition upon the operation of the synthesizer 
rests with its ability to generate the plane wave. Figure 29 shows that it 
is possible to generate the plane wave to within +j>° of phase at the indicated 
measurement position. In these curves phase is measured with respect to the 
0=0° position as shown in figure 8. Perhaps the deviation in this curve 
can be better understood by reference to figure 25. There, the array's direc- 
tivity for m = 0 indicates that the main beam has a slight asymmetry. Since 
this measurement was carried out with the 1/4” rotating microphone drive 
disassembled, reflection from that processing system can probably be ruled out. 
During these tests each driver's output was sensed to be within +0. 5° of phase 
and driven at the same output level. But these numbers were only accepted 
from the center of each speaker's output. Sufficient care was not taken at 
this time to measure accurately the entire near field pressure level and phase 
existing over the entire near field area of each driver. Variations in driver 
characteristics are believed to produce the deviation shown in figure 11a. 

Later experiments with better model drivers verified this hypothesis. Since 
all measurements at different axial locations and solid angles exhibit the 
same basic deviation from the expected wave, it is felt at this time that 
nonuniformities of each driver's near field sound distribution is the princi- 
pal basis for the majority of the discrepancy. In future experiments one 
would be wise in the purchase of commercial sound drivers of higher quality. 

As can be seen in figures 30, very good agreement between measured and 
expected phase distribution is obtained when phasing for mode 1. Apparently, 
as the solid angle ip was increased from 15° to 25° a greater divergence from 
the expected distribution occurs. With reference to the array's directivity 
for mode 1, figure 26, one can infer that the effect of driver inhomogeneities 
on the mutual impedance between drivers manifests itself in producing a more 
distorted near pressure field. This trend continues as the phase angle is 
increased between drivers to generate modes 2 and 3 as is evidenced by figures 
31 and 32 resprectively . 

It is possible to calculate the mode shapes from the first set of curves 
and measurements of the sound pressure .level at each of the 0-coordinate 
positions where phase was recorded. This essentially results in the sound 
pressure distribution as a function of 0 for an arbitrary fixed time. Once 
again the 0=0° position was chosen to correspond to maximum amplitude for 
the time selected and the other points were calculated from this reference. 

In each curve of figure 33 to 35 the expected waveform is plotted on the 
basis of the average amplitude which was recorded from all the measurements. 
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For mode 0 excellent agreement is obtained with the expected distribution, 
and as is evidenced in figure 34 the waveforms for modes 1 and 2 demonstrate 
that it is very possible to generate with the synthesizer high quality spin- 
ning waves. As may be expected from viewing figure 28, the near field inter- 
action between sound drivers induces considerable far field distorted signals 
when phasing for mode 3. In fact one should not expect a reasonably defined 
waveform for solid angles < 20°. Since this constitutes the maximum solid 
angle for measurements, data for mode 3 is not presented. In any case the 
results appear satisfactory, in the sense that at least up to mode 2 high 
quality spinning waves can be generated using only eight equally spaced 
drivers in a single ring. 


CONCLUSIONS 


This report describes research directed toward the development of a 
device capable of simulating the discrete set of eigenfunctions as produced 
by fan or compressor vehicles. The system we chose consisted of only 
eight piston-like sound drivers which permitted only the development of the 
first two spinning modes. Figures 34 and 35 clearly indicate that for 
this system the modes which were developed were of high quality, and the 
ease by which nBft and m could be changed provides an excellent means for 
parameter variation when this system is coupled to a duct to study the attenu- 
ation properties of various acoustic materials in the presence of these 
spinning waves . 

When coupling this device to a duct, it has been determined as in 
Appendix C that the most efficient location for the sound drivers would be 
in the outer 10% of the duct radius. In this position the maximum number of 
(m,y) modes could be selectively generated and propagated in the duct. 

Finally, the results of figure 26 suggest that this device could also 
be used to study the sound pressure distribution of an engine inlet by 
synthesises of its sound field. Since reliable sound measurements are diffi- 
cult in the presence of airflow at the engine inlet, the far field directivity 
patterns of the engine and synthesizer could be measured and compared, and 
the latter adjusted to give identical results on a filtered peak. 
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APPENDIX A 
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Symbols 

piston or driver radius 
amplitude of the (m,n)th mode 
duct radius 

# of rotor blades 
lo,cal sound speed 
radial phase velocity 
axial phase velocity 
circumferential phase velocity 
hub radius 

sound driver separation 
radial duct eigenfunction 
acous t i c w avenumb er 
radial duct eigenvalue 
axial duct eigenvalue 
circumferential duct eigenvalue 

harmonic of blade passage frequency (time eigenvalue) 

acoustic sound pressure 

source acoustic volume flow rate 

radial coordinate 

source ring radius 

distance from driver plane (defined in figure 8) 
source volume velocity 

# of stator vanes 
axial coordinate 
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mn 
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spiral angle in the (r,0) plane 
the Dirac delta function 
1 when m = 0, and 1/2 when m > 0 

circumferential coordinate 
radial mode index 
fluid density 
hub tip ratio (d/b) 

solid angle from origin of each driver (defined in figure 8) 
adjustable phase of each sound driver (equation 11) 

phase of (m,n)th mode 

solid angle from the origin of driver plane (defined in figure 8) 
radian frequency of compressor or sound drivers 
rotor RPM 
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APPENDIX B 


Calculation of the Allowed Order of Circumferential Modes 


The derivation for the allowed order of circumferential modes does not 
appear in any way to be unique. The procedure adopted here relies mainly 
upon a suggestion proposed by Bragg and Bridge [2] which treats the geometri- 
cal characteristics of propagating spiral waves from a ducted rotor-stator 
interaction. In figure B1 we attempt to show how the periodic redistribution 
of pressure over rotor blades due to the stator velocity defect gives rise 
to a fluctuating lift force. Now it is known (See Lowson[3]) that fluctua- 
ting lift forces gives rise to dipole sound fields and a careful observation 
of the coordinates of the wake-blade interaction zone reveals that there are 
v(# of stator vanes) interaction zones fixed at regular intervals of 2tt/v 
radians in the 0 coordinate. In each zone the dipole pressure pulse moves 
from the hub to the tip of the rotor in a cyclic fashion with a period of 
27r/vft seconds. 

As discussed in the main text, those normal modes, (m,y) for which 
equation 2 is real are propagated unattenuated in the duct. As is obvious, 
should a propagating mode be one with m > 0, then its wavevector or direction 
of propagation will be directed off-axis. Thus its lines, of constant phase 
which define the wavefront will form a helical curve of increasing pitch 
diameter along the duct axis. This condition is shown schematically in 
figure 6 and its distribution for a fixed time in figure 7. 

The spiral angle of these wavefronts formed in the (0,z) plane is deter- 
mined by the vector addition between circumferential and axial particle phase 
velocities as is expressed by equation B1 . Since all acoustic disturbances 


l 

c 





_ c e c z 

Cx = V c 0 


(Bl) 


propagate at the local sound speed, their circumferential phase velocity 
must be sufficiently supersonic to guarantee a sonic condition at the wave- 
front. Otherwise the velocity triangle in figure B2 could never be satis- 
fied. As one can see, the circumferential phase velocity is proportional 
to the rotor speed ft, number of rotor blades B, harmonic index n, and radius 
r, but inversely proportional to the circumferential mode number m. 

As is indicated in figure Bl, the spiral angle a is given by 


cx 

sm a = — 
c 0 


cx m 
27mBftr 


(B2) 
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As we have mentioned previously, we can consider the interaction of a rotor 
and stator wake to consist of v sources of pulsating pressure. By defini- 
tion, the acoustic pressure radiated from each of these sources will arrive 
in phase along all points which define the wavefront. Shown in figure 83 
is a row of stator sources radiating sound in such a manner that the* wave- 
front forms the spiral angle a in the (0-z) plane. At the point P, the 
acoustic pressure is in phase with each source by any multiple of 2tt. This, 
of course, implies that the phase relation among the sources coupled with 
varying distances between P and each source must add to be a multiple of 
2tt. For a rotor cascade consisting of B blades and rotating at Q rad/sec, 
this can be expressed as follows: 

Source Phase Rotation Phase from Increased Distance 


2 tt . nri , 2tt sin a * 0 , 

* nBft + — — *nBftr = 27rh 

c x V 

where h = 0, 1, +_ 2, . . . Equations B2 and B3 can be solved simultaneously 

for m producing the classical result first reported by Tyler and Sofrin [1] 
and given by equation B.3. 


m = nB + hV 


(B.3) 



APPENDIX C 


Effect of Source Location on Output Spectrum 


Consider the radiation of sound from a ring of simple sources with a 
specified 0 distribution into a semi-infinite rigid tube closed at z = 0 . 

This configuration is displayed in figure 2. Assuming that the sound field 
is simple harmonic, P (x,t) = P (x) e-iu)t, then the classical homogenous 
Helmholtz equation to be solved is 

v 2 P w 00 + k 2 P u (x) = o (Cl) 

with the inhomogenous boundary condition at z = 0, i.e. 

9P 

I =ikpc U(jO (r , ,0 f ) (C2) 

z=0 

where the volume velocity is specified to be Q for the particular case of 
all sources in phase and is given as 

y^O'.e') = cos (n0 1 + <j>^) 6(r-r e ) (C3) 

The remaining boundary conditions are those typically expected for a circular 
rigid tube with no reflected waves, and are 

1) 9P^/9r =0 at r = b and 3P^/9z = 0 at z = 0, / z 1 

2) P is finite as r ■> 0 and z oo 

(jl) 

3) P^ and 8P^/3 0 continuous at 6 = 0, 2 tt 

In general, the most useful technique for solving the above system is 
to make use of the appropriate Green* s function associated with the eigen- 
functions of the completely homogenous system. Essentially this offers the 
greatest flexibility should one later decide to refine the acoustic source 
function into a more realistic distribution. The most efficient method for 
determining the Green's function for this problem is by the method of eigen- 
function expansion. This technique is described particularly well by Morse 
and Feshbach [5] and is valid since the system is entirely self-adjoint. 
Applying this technique, we find that the appropriate Green's function for 
this problem which satisfied the adjoint boundary conditions, continuity at 
the source plane, and the jump condition in the first derivative is given by 


G tx,x') 
go — — J 
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ie cos(m0+(J> ) cos (m0 f +<j> *) J (k r) J (k r J 

m v nr v Y nr m my J m my 

Trt^" 1 k [J , (k b)]2 

m,y z L m+l my 


* cos k z e 
z 

ik z 

cos k z ' e Z 


ik z' 

2 0 < z < z 


z > z ' 



where 


-f 


m = 0 


£ 


m 


W 2 c 



b 2 


1) m > 0 


k = [(^) 2 - k 2 ] l/z 
z LV c' my 


k^ = set of numbers generated by — ^ 


9 J (k r) 
m v my ' 


r=b 


0 


When observation of the tube sound pressure is desired for positive axial 
displacements away from the source plane the lower z-coordinate eigen- 
function appearing in braces in C. 4 is to be used. 


At this point we can make use of Green's cookbook formula for deter- 
mining the tube's sound pressure field, i.e. 


PJx) = fff s u (x*) G u (x,x') d V' 


+ ff 


9P Cx') 

-57— - p U £ i ,) 


3G u (x,x') 

“ TP 


d S' 


Since there are no sources inside the tube's volume and the outward pointing 
normal gradient of the Green's function is zero, we only need to evaluate 
the surface integral as given in equation C.5. Inserting 

9P (x') 

P^W = // G u (x>x') — aTi — r* dr' d0' (C5) 

equations C.2, C.3, and C.4 into C.5 the tube's sound pressure field is 
found to be with the source plane at z 1 =0 


P(x,t) = P w (x) e" lt0t 


£ m k p cQ 

(2Tr)(irb 2 ) Z 
m,y 


cos(m0+(() ) J (k r) J (k r Q ) 
v Y nr m my ' m v my J 

k [J (k b)] 2 
z L m+i v my 


* 2tt 

f cos (m0 ’+<()') cos (n0 '+<(>') d0' 


This last integral is identically zero unless m = n, which implies that the 
only 0-tube mode excited corresponds exactly to the desired 0 source distri- 
bution. When m = n the value of the integral is tt and we have finally for 
the tube sound pressure 


P(x,t) 


^ k P C Q 
2irb 2 


cos (m0+ 4> m ) J m (k my r) J m (k my r 0 ) 

OCbTT 


K [J. 


m + 1 my 


(C6) 


* i(k z z - uit) 
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Figure 3 Attenuation Factor J q C^ 0 ^ r 0 *j 
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Figure 4 Attenuation Factor 







Cutoff Frequency 


Tfr O* O- 04 
04 \D OO 00 
LO O LD O 
04 tO tO Tt 




O00vDTfO4Ot0L0l^0> 
»-* o o o o OOOO 

I I I I 



23 


Figure 5 Attenuation Factor 
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Figure 6 Visualization of Spiral Wave 
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Figure 7 Predicted Waveforms in the (0,r) Plane for a Fixed Time 
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Figure 8 Coordinate System for Determining Acoustic Field of a System of n Sound Drivers 
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Figure 9 Calculated Sound Pressure Level vs. Theta 
for Mode 0 at 2000 c/s 



Figure 10 Calculated Sound Pressure vs. Theta for Mode 1 at 2000 c/s 
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Figure 11 Calculated Sound Pressure vs. Theta for Mode 2 at 2000 c/s. 
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Figure 12 Calculated Sound Pressure vs. Theta for Mode 3 at 2000 c/s. 
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Figure 13 Calculated Sound Pressure vs. Theta for Mode 4 at 2000 c/s. 



31 


Figure 14 Calculated Sound Pressure vs. Time for Mode 0 at 2000 c/s 
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Figure 15 Calculated Sound Pressure vs. Time for Mode 



33 


Figure 16 Calculated Sound Pressure vs. Time for Mode 2 at 2000 c/s. 
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Figure 17 Calculated Sound Pressure 





(tt/2) x At = 3140 rad/sec. 
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Figure 18 Calculated Sound Pre 
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Figure 19 Comparison of 1828C and Piston Driver Radiation Patterns 
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Figure 22 Assembly Drawing of 1/4" Microphone Sound Processing System 
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Figure 23 Equipment Schematic 
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Figure 24 Phase Splitting Network 
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Figure 26 Sound Driver Array Directivity for Mode m=l at 2000 c/s 
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Figure 27 Sound Driver Array Directivity for Mode m=2 at 2000 c/s 
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Figure 28 Sound Driver Array Directivity for Mode m=3 at 2000 c/s 
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MEASURED PHASE ANGLE (DEGREES) 
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ROTATING MICROPHONE POSITION (DEGREES) 


Figure 30 Measured Phase vs. Theta for Mode 1 at 2000 c/s. 
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MEASURED PHASE W/R TO ARBITRARY SIGNAL (DEGREES) 
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Figure 31 Measured Phase vs. Theta for Mode 2 at 2000 c/s 
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MEASURED PHASE W/R TO ARBITRARY SIGNAL (DEGREES) 
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Figure 32 Measured Phase vs. Theta for Mode 3 at 2000 c/s. 
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Figure 33 Measured Sound Pressure vs. Theta for Mode 0 at 2000 c/s. 
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Figure 34 Measured Sound Pressure vs. 0 for Mode 1 at 2000 c/s. 
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